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Abstract 

We perform a 3D reduction of the two-fermion homogeneous Bethe-Salpeter 
equation, by series expansion around a positive-energy instantaneous ap- 
proximation of the Bethe-Salpeter kernel, followed by another series ex- 
pansion at the 3D level in order to get a manifestly hermitian 3D poten- 
tial. It turns out that this potential does not depend on the choice of the 
starting approximation of the kernel anymore, and can be written in a 
very compact form. This result can also be obtained directly by starting 
with an approximation of the free propagator, based on integrals in the 
relative energies instead of the more usual 5— constraint. Furthermore, 
the method can be generalized to a system of N particles, consisting in 
any combination of bosons and fermions. As an example, we write the 
3D equation for systems of two or three fermions exchanging photons, in 
Feynman or Coulomb's gauge. 

PACS ll.lO.Qr Relativistic wave equations. 

PACS 11. 10. St Bound and unstable states; Bethe-Salpeter equations. 
PACS 12.20.Ds Specific calculations and limits of quantum electrodynamics. 
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1 Introduction. 

The homogeneous Bethe-Salpeter equation || is the usual tool for computing 
relativistic bound states. The first difficulty of this equation comes from the 
presence of N-1 (for N particles) unphysical degrees of freedom: the relative 
time-energy degrees of freedom. In the two-body problem, the relative energy is 
usually eliminated by replacing the free Green function by an expression com- 
bining a delta fixing the relative energy and a free 3D propagator [3- 17]. The 
exact equivalence (in what concerns the physically measurable quantities of the 
pure two-fermion problem) with the original homogeneous Bethe-Salpeter equa- 
tion can be obtained by recuperating the difference with the original free Green 
function in a series of correction terms to the 3D potential. The dependence 
in the relative energy of the 4D Bethe-Salpeter equation is thus transformed 
into a series of higher-order correction terms to a 3D equation. This series of 
'modified ladder terms" is to be combined with the series of "crossed terms" 
which gives the Bethe-Salpeter kernel as a sum of terms corresponding to ir- 
reducible graphs. We do not know if the result converges, and, if yes, which 
will be the level of approximation obtained with a given truncation. This will 
of course depend on the specific problem studied. The mutual cancellation be- 
tween the higher-order " modified ladder terms" (generated by the 3D reduction) 
and the "crossed terms", which occurs at the one- body hmit (l|, |l|, |l|l IS an 
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encouraging fact. The starting homogeneous Bethe-Salpeter equation (for the 
Bethe-Salpeter amphtude) is deduced from the inhomogeneous one (for the full 
propagator) by postulating the presence of a bound state pole, and is thus a 
priori valid only for bound states. However, the scattering amplitude built with 
the 3D free propagator and potential turns out to be correct on the mass shell. 

This kind of 3D reduction of the two-fcrmion homogeneous Bethe-Salpeter 
equation has been performed by many authors [3-17]. The resulting equations 
are theoretically equivalent at the limit of all correction terms included. Of 
course, in practice, the series generated by the 3D reduction and the series 
giving the Bethe-Salpeter kernel itself will both have to be truncated and the 
estimation of the level of approximation obtained is not straightforward. 

A less often used method is based on the replacement of the Bethe-Salpeter 
kernel by an "instantaneous" (i.e. independent of the relative energy) approx- 
imation. In this case, the resulting 3D potential is not manifestly hermitian 
(for each fixed value of the total energy on which its depends). In Phillips and 
Wallace's method, the starting approximation is tuned, order by order, is such 
a way that it becomes the final 3D potential . In a recent work, we also got 
an hermitian 3D potential by performing a supplementary series expansion at 
the 3D level and combining it with the first 3D reducing expansion [ p^ . 

In the three-fermion problem, a new difficulty comes from the unconnect- 
edness of the two-body terms of the Bethe-Salpeter kernel, which are in fact 
the more important terms and often the only ones to be considered. In a con- 
straining propagator-based reduction, the " approached" propagator would have 
to be a 3D propagator with two constraints (fixing the two relative energies). 
However, when an unconnected two-body term lies between two of these ap- 
proached propagators, we get a constraint too much. In Gross' spectator model 
[ pO| the transition matrix element T is approached by putting on their 
positive-energy mass shell the fermion lines which are not going to be inter- 
acting before and after. This means that each three-fermion free propagator 
will be " approached" in a way depending on its position in the expansion of T. 
This leads to a set of Faddeev-like equations, which is not equivalent to a single 
equation with a potential, unlike in the nonrelativistic case. 

Recently, we succeeded in adaptating our two- fermion kernel-based reduction 
to the three-fermion case. The 3D potential was a complicated sum of terms 
of various origins, although the first-order approximation was still manageable 
When working on an application of this 3D reduction method, we revisited 
our two-fermion kernel-based 3D reduction, and found an alternative to the 
second series expansion (the one performed at the 3D level to get an hermitian 
potential) . After combining this series expansion with the first one (given by the 
3D reduction itself), we found that the starting instantaneous approximation of 
the Bethe-Salpeter kernel disappears from the final 3D potential. The result is 
in fact a compact expression of the potential that Phillips and Wallace compute 
order by order [ |T8| . Furthermore, it can more directly be obtained by a new 
integrating propagator-based reduction method, in which the relative energy 
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is integrated on, instead of being fixed by a 5— fonction (or constraint). Here 
again, the physical scattering amphtude computed at the 3D level is also correct. 
In fact, the 3D full propagator turns out to be the retarded part of the 4D full 
propagator with the initial and final relative times put to zero, as proposed by 
Logunov and Tavkhelidze . 

This integrating propagator-based reduction (and the Logunov and Tavkhe- 
lidze's reduction can easily be adapted to the three- fermion problem, 

as the unconnectedness of the two-body parts of the Bethe-Salpeter kernel is 
not a difficulty anymore when integrals on the relative energies are used instead 
of (5— fonctions. Furthermore, the method can be generalized to a system of N 
particles, consisting in any mixing of bosons and fermions. The 3D free prop- 
agator is given by the integral of the product of the N free Green fonctions 
with respect to the N-1 relative energies and consists in two terms: a term in 
which all particles have positive free energies, and a term in which they have 
all negative free energies (a slightly different 3D reduction, with simpler first 
terms, could be obtained by removing this part of the 3D propagator. We shall 
consider both options below). The absence of mixed free-energy signs preserves 
the 3D equation from continuum dissolution ||2^, ^ |2^, The 3D scatter- 
ing amplitudes are again correct when all particles are on their positive-energy 
mass shell (which has no measurable consequence when iV > 2 ) , but also in more 
complicated processes, like the scattering of a free particle by a bound state. 

There are many ways of performing a constraining propagator-based 3D 
reduction of the two-body Bethe-Salpeter equation, according to the choice of 
the constraint and of the 3D free propagator. Our integrating propagator- 
based reduction approach, in contrast, leads to unexpectedly unifying results: 
independently of the choice of the initial instantaneous positive-energy kernel, 
the final 3D equation is always the same. Furthermore, it can also be obtained 
by Phillips and Wallace's method and by the quite different Logunov and 
Tavkhelidze's method while our integrating propagator-based reduction 

provides a more compact form for the 3D potential by recombining the series. 
We do not pretend that this 3D reduction converges more rapidly that the more 
usual constraining propagator-based reductions (we try a comparison below). 
The important result consists in the fact that it can easily be generalized to 
more than two particles, giving a 3D equation for a wave function with a 3D 
potential, which must necessarily be truncated but can be indefinitely improved. 
By comparison. Gross' spectator model for three particles |20[ |l| (which can 
be considered as a constraining propagator-based reduction) leads to a set of 
coupled equations for partial transition operators, which can not be reduced to 
a single equation. 

As a first test of our reduction method, we choose the problems of two and 
three fermions in QED, using Fcynman's gauge. We compare with a constraining 
propagator-based reduction method in the two-fermion case, and with Gross' 
spectator model in the three-fermion case. In actual calculations. Coulomb's 
gauge would be a better choice than Feynman's gauge, but it is less suited to a 
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presentation of our method, as the corrections to the instantaneous approxima- 
tion and the three-body terms are then of higher order, competing with a lot of 
other contributions. 

Two natural requirements in the writing of phenomenological equations for 
systems of mutually interacting relativistic particles are Lorentz invariance and 
cluster separability: if we split the system into clusters by "switching ofP' selec- 
tively the interactions between the clusters or by by letting the distance between 
them become infinite (this is not quite the same thing), we want to get inde- 
pendent sets of correct covariant equations. It is not easy to satisfy Lorentz 
invariance and cluster separability together, and, moreover, to include contin- 
uum dissolution-preserving operators. In our 3D reductions of Bethe-Salpeter 
equations the approach is different: the original inhomogeneous Bethe-Salpeter 
equation (from which the homogeneous one is derived) is manifestly Lorentz co- 
variant and cluster separable. The 3D reduction is performed in an unspecified 
(at first) reference frame and the individual terms of the series giving the 3D 
potential are not covariant, some of them being not cluster separable either. If 
we choose to work in the global center of mass frame, all terms can be made 
formally covariant, but this reference to the total momentum spoils the cluster 
separability. However, if we see the 3D reduction only as a tool to calculate 
the bound state spectrum, the exact results can in principle be indefinitely ap- 
proached by including more and more higher-order contributions. Moreover, it 
turns out that this same 3D equation is still valid for the various physical scat- 
tering amplitudes. The computation of electromagnetic form factors can also 
be done by computing electromagnetic scatterings. What is then the practical 
meaning of the cluster separability requirement? the selective "switching off' 
is clearly a purely mathematical operation which does not correspond to some- 
thing physical, while the problem of clusters going to infinity must be treated 
as a scattering problem. A practical aspect of the Lorentz covariance / clus- 
ter separability requirement could be the use of, for example, two-body bound 
state or scattering data in a three-body context. When used directly, the two- 
body data must always be extrapolated off the mass shell, in a model-dependent 
way. In our 3D reductions, we could use two-body data indirectly, through the 
adjustment of coupling parameters. 

In section 2, we present the usual constraining propagator-based reduction of 
the two-fermion Bethe-Salpeter equation, followed by our kernel-based reduc- 
tion transformed into an integrating propagator-based reduction, and we com- 
pare the first terms of the 3D potentials obtained by both methods. In section 3 
we generalize this integrating propagator-based reduction to N-particle systems 
(fermions, bosons or mixings of both). In section 4 we apply our method to sys- 
tems of two or three fermions exchanging photons, in Feynman's and Coulomb's 
gauges. 
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2 The two-fermion problem. 

2.1 Constraining propagator-based reduction. 

We shall write the Bethe-Salpeter equation for the bound states of two fermions 
i as 

$ = G'"X'$, (1) 

where $ is the Bethe-Salpeter amphtude, function of the positions xi,X2 or of 
the momenta pi , p2 of the fermions, according to the representation chosen. The 
operator K' is the Bethe-Salpeter kernel, given as a factor of the kernel of an 
integral equation in momentum space by the sum of the irreducible two-fermion 
Feynman graphs. The operator G'° is the product G'iG'2 of the two individual 
fermion's dressed propagators 

G? = ^ ^ ^. (2) 

where rrii is the mass of the fermion i and (p^) the renormalized self-energy 
function, which becomes proportional to (7^ • pi — m^)^ near the mass shell. 
We shall need a Bethe-Salpeter equation in terms of the free propagator 

G° = G?G°, G° = — ^ (3) 

■^i ■ Pi - nil + «e Pio - ni + lehi 

where the are the Dirac free hamiltonians 

hi = di .pi + Pi rrii (i = l,2). (4) 
In order to transfer the self-energy part of G'" to the kernel, we shall write 

K = K' + (5) 

with 

S = (G")-i - (G'")-i = Si(G^)-i + S2(G;)-i - S1E2 (6) 
and get the Bethe-Salpeter equation in the form 

$ = G°if $. (7) 

We shall define the total (or external, CM, global) and relative (or internal) 
variables: 

X^-{XI+X2), P=Pl+P2, (8) 

X = Xi-X2, P=^(Pl-P2)- (9) 

and also 

S = hi + h2, E = Ei+E2, = a/^= (p?+m2)i (10) 
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We do not specify the reference frame in which we write noncovariant quantities 
Uke hi or Ei. Our 3D reduction will in fact be frame-dependent. Practically, 
we shall choose the global rest frame of the 2- (in this section), 3-, or N-particle 
system. In the N-particle problem we shall not try to boost the two-body 
kernels between the global rest frame and the (virtual) two-particle rest frames: 
these relativistic effects will be taken into account (if desired) by the inclusion 
of the higher-order terms of the series generated by the 3D reduction, since 
this series, when untruncated, leads to the same measurable quantities as the 
starting covariant Bethe-Salpetcr equation. 

If we consider the contributions of the poles of 

= 1 1 Pip2 (11) 

|Po +Po - hi+ iehi ^Pq - po - h2 + ieh2 

in an expression like KG'^K, we must perform an integration with respect to 
Po- If K is instantaneous, we get 

/ dpoG°(po) = -2i7rTg°(3,P2, 9° = -5 . „ (12) 

J Po - S + lePo 

where ^ ^ ^ 

^ = 7;in + 7-2), Ti = = (13) 

can also be written 

r = A++-A--, A^^=A1A^2, = (1^) 

When K is not instantaneous, we must add the contributions of its singularities. 
Furthermore, in the residues of the poles of G° we must take K at pio = h\ 
(i.e. at +Ei in the subspace built with the h\ = +Ei eigenstates of hi , and at 
—El in the complementary subspace) or at P20 = ^2j according to the chosen 
integration path and to the sign of r. We shall perform a 3D reduction based 
on the replacement of the free propagator G° by a carefully chosen expression 

G^(po) = -2nrTS{pio-hi)g" f]if52 (15) 
combining the 3D propagator 5" with the constraint 

S{pio-hi) = At6{pio-Ei) + A^S{pio+Ei) = d{po-si), si = y + Zii- (16) 

The operator r has a clear meaning in the basis built with the free solutions: it 
is -hi for hi,h2 > 0, -1 for /ii, /12 < and zero when they have opposite signs. 
It comes from the dependence of the po integral on the signs of the iehi. 

This choice of has three merits: a) It leads directly to a simple equa- 
tion (Salpeter's equation without higher-order correction terms to the potential) 



7 



with an instantaneous approximation of the kernel. The higher-order correction 
terms (necessary to get the correct one-body limits) will then be the consequence 
of the dependence of the kernel in the relative energy, b) In the two-fermion 
plus potential problem and in the three-fermion problem, the operator r pre- 
vents the mixing of asymptotically free fermions with opposite energy signs, 
which is the origin of the continuum dissolution problem [^5[ |2^, ^ c) It 
preserves a particle-antiparticle symmetry, which is a characteristic feature of 
relativistic theories. We could also replace r by A++, as the A part does 
not contribute much in practice. One of the fermions would then be on its 
positive-energy mass shell, as in Gross' spectator model There is an infin- 
ity of other possible choices, like fixing pq at its on mass shell value (/ii— /i2)/2 , 
instead of fixing piQ at its on mass shell value hi , treating thus both fermions 
in a symmetrical way. With Sazdjian's propagator one gets an explicitly 
covariant equation (we recently combined this propagator with a kind of co- 
variant positive-energy projector [0). With Lepage ||], we get directly a 3D 
equation written in Schrodinger's form. An important feature is the presence or 
absence of a continuum dissolution-preserving operator like r or A++. Without 
such an operator, the propagator could not be used beyond the two-body prob- 
lem, including the two-body in an external potential problem. Our constraining 
operator contains thus the operator r (we shall also examine the A++ choice) . 
Besides that, the principal reason of our choice (|l^) is to allow for an easy term 
by term comparison with our integrating propagator-based reduction below. 

Let us now write the free propagator as the sum of the zero-order propagator, 
plus a remainder: 

G" = + G". (17) 
The Bethe-Salpeter equation becomes then the inhomogeneous equation 

$ = G^K^ = {G^ + G")K^ = * + G^A'$, (18) 

with 

* = G^K<^ (= G\G^)-^^). (19) 

Solving (formally) the inhomogeneous equation ( [l^ ) with respect to $ and 
putting the result into ([l9|), we get 

* = G'^A:(1 - G^KY^^ = G^K'^^ (20) 

where 

K'^ = K{1 - G^K)-^ = K + KG'^K + ... = (1 - KG^Y^K (21) 

obeys 

K'^ =^K + KG^K^ ^ K + K'^G^K. (22) 

The reduction series ( ^l| ) re-introduces in fact the reducible graphs into the 
Bethc-Salpeter kernel, but with G^ replaced by G^. Equation ( pO| ) is a 3D 
equivalent of the Bethe-Salpeter equation. 
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The relative energy dependence of eq. (EG) can be easily eliminated: 



* = ^(po-si)V (23) 

and tp obeys: 

^=-2i7rg°T J dp'odpo6ip'o-si)pip2K^ip'o,Po)Sipo-si)i^. (24) 
Using the identity ■!/) = t^'0, we can write 

iP^g^TV^jj, V = -2iTTT^ f3l(32K^{si,Sl)T^, (25) 

/?i/32i^^(si,si) = / dp'odpo6ip'o-Si)PiP2K^{p'o,Po)Sipo~s,). (26) 



Note that we write {p'q,Po) but (si,si), as we keep si in operator form. This 
operator can be diagonalized in the spatial momentum space by using the A'^ 
projectors. The eigenvalue will depend on the position of si in the formula: 
the eigenvalue of the first si in ( p6[ ) will be built with the final momenta and 
that of the last si will be built with the initial momenta. 
The inversion of the reduction is given by 

$ = (l-G^K)-^-^ = (1 + G^X'^)* = {l + G^K^-G^K^)-^ = G^K'^^ 

(27) 

or, explicitating the relative energy: 

^{p'^) = G\p'^)K^{p'^,s^)^. (28) 

The splitting of G^ into two terms containing a 5 is the origin of unphysical 
singularities in the terms of when the argument of the delta vanishes on 
the singularities of K. When the full is computed, the singularities of the 
different terms cancel mutually. When is truncated, some of the unphysical 
singularities have to be removed by hand |20| . 

This 3D reduction can also be described in terms of transition operators. The 
4D transition operator corresponding to our modified Bethe-Salpeter equation 

©is 

T = K + KG"K + --- (29) 

and can be obtained by keeping only the G^ part of G° in it. We have 
also 

T = Kil-G^K)-'^ = K{1-G"K-G^K)-^ = K{l-G"K)'^il-G^K{l-G'^K)-^)-^ 

= K^il-G^K^y^ ^K^ + K^G^K^ + ■■■ (30) 
so that the 3D transition operator 

y + y g° r F + • • • (31) 
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is also given by 

= -2nrT^ PiP2T{si,si)t\ (32) 

The operator T and the operator T' corresponding to the original Bethe- 
Salpeter equation (|^) are related by the common full propagator G : 

G° + G°TG" = G = G'" + G'"T'G'° (33) 

with 

G'° ^ G"{1-^G° )-i = G? (1 - Si G? )-iG^ (1 - S2 G^ )-i. (34) 

We see that the 3D transition operator is a constrained form of that of field 
theory. When both fermions are on their positive-energy mass sheUs, the op- 
erators T{si, si),T,T' , become equal and T^^ becomes proportional to the 
physical scattering amplitude. This was not guaranteed a priori, as our original 
two-fermion Bethe-Salpeter equation (0) was valid only for bound states. Our 
3D equation (|2^) is a bound state equation too. To include the scattering states 
we should add an inhomogeneous term, or write the equation in the form 

{Po-S)4> =tV4'- (35) 
2.2 Kernel-based reduction. 

A few operators in this kernel-based reduction will be given the same name as 
similar operators used in the constraining propagator-based reduction, although 
being not identical to them. In order to avoid any confusion, we shall overwrite 
these operators with a ~ . 

If the Bethe-Salpeter kernel were instantaneous (i.e. independent of pq), 
one would get Salpeter's equation by integration with respect to po In 
the realistic case of a non-instantaneous kernel, it is possible to compute the 
bound state energies at the 4D level by perturbations around an instantaneous 
approximation of the kernel CO, 31, 1^. Here, we want to build a 3D reduction 



around an approximation K of the Bethe-Salpeter kernel. Let us write 

K = K° + K'^. (36) 
The Bethe-Salpeter equation becomes 

$ = G°i^°$ + G"i\:^$ (37) 
<^ ^ {l-G°K^)-^G°K°<^ ^ G^K^<^ (38) 

with 

G^ - G° + G°if ^G° + • • • = G" + G^«. (39) 
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If we now specialize K'^ to an instantaneous positive-energy kernel [K'^ = 
A++/3i/32-f^°A++ and is independent of po), we get 

= (40) 

with 

= A++ y dpo $(po) , V" = -2z7r/3i/32 K° , (41) 

g^'' = A++ 1 dp;,dpo G'^^(pE„Po) /3i/32 A++. (42) 

The interaction term (5")^^( g° + g^^ ) is not hermitian (for a fixed value 
of the total energy on which its depends), so that we do not know if equation 
( ^0| ) will have a real energy spectrum. Let us write this equation in the form 

- (1 + .gOT^«) , T^"" = (.9°)"\9'^''(g°)"' (43) 

in order to perform the transformations 

(1 + 5"r^^)- V = (44) 

= + 1 - (1 + .gOT^^)-i ] = gO [ + r^^(l + .9"r^^)-i ] 

= 5° [ < ° > + < i4:^(l - G^K^)-^ >{l+g° < K^{1 - G^K^)-^ >)-^ ] 

(45) 

with the definition 

< A>=-iA++(.gO)-i / dp^dpoGO(p^)A(p;„po)G"(po)/3i/32A++(.gO)-i. 

(46) 

Let us expand g'^ into a 4D operator: 

= > 5° < (47) 

defined as (writing the relative energy arguments, the other momentum argu- 
ments P and p remaining local): 

A++ 

G'(p'^,Po) = G\p'^) /3i/32 — G"(po) (48) 

— ZITT g 



and write ( 45 ) in the form 

= < ifO + x«(l - G'°i^^)-i(l + G'^ii:^(l - G^i^^)-!)-! > (t> 

^ g" < K°+K^il-G°K''+G'K'')-^ > <j> = g" < K°+K"{1-&K'')-^ > , 

(49) 

with ^ 

G«(p^,Po) = G°(po)'5(p;, -po) - G'ip'^po). (50) 
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It is easy to see that AG^ — G^A — whenever A is an instantaneous positive- 
energy operator hke K^. We can then reconibine K'^ + = K in the 3D 
equation, which becomes 

(t)^g'^V4>, V^<K'^>, ^ K{1-G'^K)-\ (51) 

Computing the scattering matrix element as in the preceding subsection gives 

f^^ = V{1~ g"vy^ = <T > . (52) 

On the positive-energy mass shells Po—E[+E2 = Ei+E2 the operators (5°)"^ 
of ( ^6| ) vanish and we remain with the residues of the corresponding poles in 
the integrations with respect to p'^ and po, so that 

<T>= -2i7rA++pip2 T{p[o= E,, p,o = E, ) A++ (53) 

and we have again the correct physical scattering amplitude ( T is written here 
with operators as arguments, as in (p6|)). 

Equation ( ^2| ) can be inverted to give V in terms of < T >: 

V ^<T > {l+g° <T >)-^ =< > . (54) 

In fact, it is possible to use equation ( p2[ ) to define the 3D reduction. The 
corresponding 3D full propagator 

g = 5"A+++.g°f3V (55) 

is then proportional to the integral of G with respect to the initial and final 
relative energies, between two A++ projectors. This manipulation preserves 
the position of the bound state poles and the physical scattering amplitudes. 
In configuration space, it means that we take the retarded part of G at equal 
times. The transformation of the inhomogeneous Bethe-Salpeter equation into 
an equation for g was the starting point of the 3D reduction of Logunov and 
Tavkhehdze 

Phillips and Wallace's 3D reduction is a kernel-based reduction in the — 1 
subspace. Their way of symmetrizing the interaction term consists in tuning 
in order to make g^^ vanish, considering thus g^^ = as an equation in if" 
to be solved order by order. If we do that here in the A++ ~ 1 subspace, we 
see that it implies that the operator at the left of (H) is the identical operator. 



The comparison with (51) leads to 



-2iTT(3i(32K° =< K{1-G^K)-^ > . (56) 

Our 3D equation ( ^l| ) is quite general, as it does not depend on the initial choice 
of K'^ anymore (at the 3D level, it is of course still possible to perform a pertur- 
bation calculation starting with < A'" >, using (49)). This 3D reduction is also 
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identical to that of Phillips and Wallace , obtained by tuning an unspecified 
instantaneous kernel, and to that of Logunov and Tavkhelidze , obtained by 
taking the retarded part of the full propagator G at equal times. Our own 3D 
reduction gives a compact expression of the 3D potential, by recombining the 
series. 

The independence of our 3D equation on the choice of K'~' suggests that it 
could also be obtained by a new kind of propagator-based reduction, in which 
the relative energies are integrated on instead of being "constrained". We shall 
now present this new approach, and we shall see in section 3 that is can be 
easily generalized to N-body systems. 



2.3 Integrating propagator-based reduction. 

We shall now perform directly an integrating propagator-based reduction, in- 
spired by our kernel-based reduction above, but we shall work in the larger 
= 1 subspace, defining 

G'{p'o,po) = G'^ip'o) Pih G°{p,) (57) 

—Znrg" 

G° = G^ + (58) 

$ = G°K<^ = {G^ + G^)K<^ = $ G^if$, (59) 

with ^ _ _ _ _ 

* = G^K'^ = G'K{1 - G^K)-^ * = G^K'^^ (60) 

where 

K'^ = K{1 - &K)-^ = K + K&K + ... = (1 - K&y^K. (61) 



When we explicitate the relative energies in equation (60), we get 



*(Po) - ^°(^'o)/?i/?2 3^^ j dp',dpoG'>{p'o)K{p'„po)Hpo) 
= G'{p'^)fi,f5,—^ I dp'oHp'o) 



G°{p'^) Pif32 / dp'^dpo G\p'^)K^{p'^,po) ^{po)- (62) 



Defining now 



= dpo$(po) (63) 
<A>^ ^r\g")-^J dp'odpoG°{p',)Aip'o,po)G'{po)Pip2T^ {g")-' (64) 
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so that 

= > rg° < (65) 
and writing (|6^) in terms of 0, we get 

= .g" r < k'^ > (j). (66) 

If we perform again our kernel-based reduction of subsection 2.2, but now in 
the T^ — l subspace and using the definitions ( |63[ - |65| ), we get 

X^g°T[<K°>-T<K°>T + T<K'^>T]x, X^T(j) (67) 

and A'" will disappear if it does not connect the (++) and the (- -) subspaces, 
leading to our equation (p6|). Phillips and Wallace's potential, which must 
reappear after the symmetrizing transformation, will be < K'^ >=< > . 

If we prefer to work in the = 1 subspace, we can recover the results of 
subsection 2.2 by simply replacing r by 

2.4 Comparison of the lowest-order terms. 

With the integrating propagator-based reduction, the final 3D potential will be 
V ^ <K> + < KG"K > + ■■■ 
= < K > + < KG^K > - < K > Tg° < K > +■■■ (68) 

At order 4 in the coupling constant, we should keep the ladder term and the 
first crossed term in < K > and only the ladder term in < K G^K > . The 
contributions of the higher-order terms should a priori decrease with the number 
of G^ : once the leading contributions containing the poles in Pq — {Ei + E2) 
have been removed, we remain with the smaller contributions coming from the 
residues of the poles of K or from its (-1--) and (--I-) components. 

For comparison with the usual constraining propagator-based reduction of 
subsection 2.1, we shaU spht into G^ +G" in the definition of <A> 
(it must be noted that this splitting introduces unphysical singularities [|l5[ ^ ) . 
Keeping only the terms with zero or one G^, we get 

V ^<K + < KG"K 

+ [1-<K rg"] <K>'^^ + <K [l - rg" < K ] (69) 

where the index R or S indicates which part of G*' has been kept. The two 
first terms are also the two first terms of the constraining propagator-based 
reduction. The two other terms would not contribute at order 4 in a pertur- 
bation calculation beginning with < K . We see the kind of rearrange- 
ments by which the constraining propagator-based reduction and the integrat- 
ing propagator-based reduction will finally lead to the same energy spectrum. 
At order 2 only in the coupling constant, we would get 

V fii < K + < K + < K (70) 
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and checking if the addition of the two last terms brings ( |70| ) nearer to ( |69| ) 
could help to choose between both approaches in a specific problem. 

The 3D free propagator is rg^ in the constraining propagator^ as in the 
integrating propagator-based reductions. The 3D potentials V and V are both 
hermitian and both equations are equivalent. We must however not conclude 
that these potentials are equal: the dependence of the potential on the total 
energy allows for an infinity of equivalent equations. 



3 N-body problem. 

3.1 Generalization of the integrating propagator-based re- 
duction. 

The integrating propagator-based reduction is a good candidate to a double 
generalization: from two to N fermions and/or from N fermions to any system 
of f=0,l,...N fermions with b = N-f bosons. The propagators of the fermion i 
and the boson j are respectively 

G° T— — r^^' (71) 

Pio - hi + lehi 

with (jj =±1. The free propagator for a system of /=:0, l...iV fermions 
and b = N-f bosons will be 

Lj — Lt]^ ... Ljj Lry^]^ ... Lj^ \''^) 

and the corresponding 3D propagator will be proportional to 

dpnG"{po) = / S {Pq - y^^pio)dpio...dpNo G°(pio, ■■■Pm)- (74) 
•' i=i 

We shall perform the integral in pm by replacing it by {Po — ^f^2 Pio) i^i the 
first propagator. Each other piQ will then appear twice: in and in the 
corresponding G^, and the integral with respect to pio will be zero unless hi 
(or (Ji ) and hi (or ai ) have the same sign. The final result will be 

dpoG^'ipo) = tl^^^rg'p (75) 



with 



•9° = p ^—-17 ' S = EiA+-A-). p = Pi...Pf , (76) 
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A±-Af...A±, r = A+ + (-)^'+iA- , (77) 

N 

E = Y,E^, Lu ^ 2''Ef+i...EN (78) 

for / 7^ 0, so that 

= p ^t^. + i-y^' p . ■ (79) 

il) — + Fq + E — le 

When / = 0, (bosons only), we have no r and no /3 (we can replace them by 
1 in (||)) and 

_ 1 1 _ 2E 

^ ^ Po-E+ it ~ Pq + E -ie ^ P2 ^ E^ + ie ' ^ 

The Bethe-Salpeter equation for N particles can still be written 

$ = G"A'$ (81) 

with 

K ^ K' + E = (G")-^ - {G'")-^ (82) 

where K' will be given by a combination of irreducible 2 < n < A^— body 
irreducible kernels. For N=3: 

K' = K[, (G^O)-i + K'^^ {Crr' + K',, (G^°)-i + K[,^. (83) 

For A^ > 4 the writing of K' becomes more complicated, because of the pres- 
ence of commutating kernels like and Ki^^ which would lead, without cor- 
rections, to an overcounting of some graphs in the expansion of G. We examine 
this problem elsewhere |Q . 

Let us now define a N-body operator G^ : 

which is such that 

dpoG^(p^,Po) = G\p'^)Pt^[3 = T^G\p'^), 



I- 



dp'oG' {p'o,Po) = t'G"{po). (85) 
Performing a 3D reduction as in subsection 2.3, we get again 

(j) ^ g^T <K^ > (j) (86) 
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with the definitions 

<j> = T^/Zi j dpo'i'ipo) (87) 
= K{1 - G^K)-\ & = -G^ (88) 

and 

<A>^ (Z^^^)lv3T ^ / dp',dpoG'ip',)Aip',,po)G'{po)f3^ (89) 
so that 

= > rg" < . (90) 

As in the two-fermion problem, we can choose to work in the = 1 subspace 
instead of the r^ = l subspace, by simply replacing r by A+. In this case also 
it is possible to start the 3D reduction by taking the retarded part of the full 
propagator at equal times p4[ . 



3.2 Scattering. 

As in the two-fermion case, we have T^^ =< T > . This implies that T^^ 
is equal to the physical scattering amplitude when all particles are on their 
positive-energy mass shells in their initial and final states. In actual scattering 
experiments, however, the initial state consists in two clusters only. 

As shown already in |2^, the 3D equation remains valid in these cases. The 
scattering of a bound state of fermions (23) by fermion 1 to a bound state of 
fermions (12) plus fermion 3, for example, can be described by a part of the full 
propagator G : 

G = G3 G12 Ii2,23 G23 Gi + ■ ■ • (91) 

where T12.23 contains neither initial (23) nor final (12) interaction, these inter- 
actions being included in G23 and G12 respectively. The full (23) propagator 
can be expanded as Q 

G23(^'23,P23'P23) = 3'23 (^^23 , P23) ^ ^ —^23(^23,^23) + ' ' ' (92) 



where we isolated the pole at P230 — E23 — y P23 + P23 , defining $23 as $^3 /32/33 . 
The full (12) propagator can be expanded similarly. The N-fermion 3D full prop- 
agator is, according to (|8S|): 

9 = ^_2il)N-i / dp'odpo G{pIpo)PA+. (93) 

Taking TV = 3 and putting (|9^) into (^3|), we see that the integration of 
^'23(-P23,P23)Aj3 with rcspcct to P230 wiU give 023(-P23,P23) while in the in- 
tegration with respect to pio we shall isolate the residue of the pole of G" at 
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pio = Ei (and similarly for the final state). This leads to 



9 = -<Pl2^ ^ ^ T? , ■ [*12Tl2,23$23]3i ^ ^ , ■ '/'23 

-Po - £-12 - Ez + ie^ '-^^ Pa - £-23 - Ei+it 

(94) 

where 

[$12^12.23^23] 31 = j dp'^dpo 5{p'^^~ E^)Ti2M, Po)5{j>w~ E^) . (95) 

If we now compute g at the 3D level, using the free propagator K^g^ and the 
potential V, we get 

<? = A3+gi2ri32%523A+ + ••• (96) 
while ^23 is obtained from G23 through (^: 



;?23 = ^ 023 F ^ F , '^'^3 + • ■ • (97) 

^TT — £/23 — ii/i + »e 



and similarly for 312. If we identify ( |94| ) to (^6|) in momentum space, we get 
the following relation between the physical scattering amplitudes computed at 
the 3D and at the 4D levels: 

<l^t2 T?2% = {-2tn f [¥12 Ti2,23 $23] 3^ /3i (98) 

the momenta being such that Po = £'i2 + £'3 = £'23+£'i- 

Here, $23 depends on P23,P23 while (/)23 depends on P23iP23- We could 
compute 023 at P23 = 0, then compute $23 from it, boost $23 to P23 7^ 0, 
and use it in the right-hand side of (p8|), or we could try to compute directly 
023 at -P23 7^0, and use it in the left-hand side of (|9^). 

We shall stop here these considerations about the scattering, which is not 
the main subject of this work and would be better studied on a specific example. 
We only wanted to show how the 3D equation not only gives the bound state 
spectrum, but is also valid for the scattering amplitudes. 

3.3 Reduction of the Dirac spinors. 

Let us first consider a fermion in a positive-energy potential: 

(po-h)^ ^V'i', V = A+VA+. (99) 
Using the A+ projector: 

(po-£^)*^ = *+=A+*. (100) 
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We shall apply a virtual boost on 



^0 ^ l + (7-^)(7-") ^+ 



with 



n=-(i;,p), M=(l, 0). 
m 



In terms of 5*", equation (100) becomes 



1 + (7 • u ) (7 • n) 1 + (7 • ") (7 • M ) ,^ ,0 



(101) 



(102) 



V 



^J{n + n)2 



Writing 



we get 



1 + /? / 2 m 



2£: 



£^ + m y/2m{E + m) 



^0 = 



2E 



E + m \ E + m 



2E 



(103) 



(104) 



(105) 



where v is the large-large part of V. Similarly, for f fermions and a positive- 
energy integrating propagator-based reduction, we have 



iPo-E)ip 



f 

n 



2E, 

Ei + m 



f 

n 



2E, 
Ei + m 



ip. 



(106) 



4 Example: two and three fermions in QED. 

Our aim in this section is only to show how the integrating propagator-based 3D 
reduction method works for a system of > 2 particles, since, as in any actual 
problem, specific complications are encountered (choice of the gauge, estimation 
of the order of the contributions for the exchange of zero-mass quanta, detection 
of the cancellations...). In Coulomb's gauge the relative energy dependence of 
the Bethe-Salpeter kernel lies in the transverse part, the contributions of which 
are typically smaller by a factor . This makes the calculation of the first-order 
contributions (up to a^) more easy, but does not provide a good illustration 
of our 3D reduction method, the three-body contributions to the 3D potential 
being then of order > a'^ and competing with a lot of other effects. We shall 
therefore work in Feynman's gauge until subsection 4.4. 
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4.1 Two fermions. 



For two fermions of charges Zi, Z2 the Bethe-Salpeter kernel is 



2iZiZ2 1 
(27r)3 fc2 + ie 



(71-72) 



2i Z1Z2 



n 2 



(27r)3 



dk 



fc'2 



fc2 



7r72^ G?(<zi)G°(g2^) 71^72, 



(107) 



We wrote explicitly the ladder term and the first crossed term K'~^ . The 
meaning of the variables can be read on figure 1, which represents 
<K^>,<K^ > and < K^G^K^ > = < K^G°K^ > - < > g° < > . 
The results of subsections 2.2 to 2.4 will be used, and specialized to the kernel 
(f07). If we choose to perform the 3D reduction in the A++ = 1 subspace, we 



> {p',p) 



2iZiZ2 
72^ 



A++(p') (1 - di-a2) A++(p ) / ip\p) (108) 



w., .^ {Po-E[-E^){Po-E,-E2, , , 



-2i7r 



k'^ + ie pio-Ei + ie P20-E2 + ie 



{E[-EiY -k^ + ie 



2\k\ 



R 



(109) 



^ ^ {E[-E,~\k\){P^-E^-E2) _ {E[-E, + \k\){Po-E[-E'^ 
Po-E[-E2-\k\+ie Po-Ei-E!,-\k\+ie 

where the first part of / corresponds to < and contributes in to 

the energy. The two terms of R correspond to < and < 

respectively. At lowest order 



R « iPa-Ei-E2) + iPo-E[~E!,). 



(Ill) 



By sheer power counting (pc), R should contribute in a'^ ( |^|~^ brings a factor 

and the difference of the energies a factor , to be combined with the 
of Coulomb's potential). However, this dominant contribution, if isolated, is in 
\k\~^ and makes the future integrations with respect to k diverge at origin. A 
more careful estimation is thus necessary, and shows that R contributes in fact 
in a^log{a). This divergence of the dominant contribution is a consequence of 
the zero mass of the photon. From now on, we shall estimate only the (pc) order 
of the contributions, keeping in mind that the real order could be lower. 
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The expression of / is apparently asymmetric for the permutation of the 
fermions, but it can be shown that it is in fact equal to the expression obtained 
after permutation or the fermions, on by directly closing the integration paths 
counterclockwise. The differences between [— [{E[ — Ei)^ — k^]~^ and 
[ (i?2 — -£^2)^ — are in (pc) (these three denominators are equal in the 
equal mass case if we work in the two-fermion rest frame) . 

A third way of computing / consists in writing the photon's pole as 



1 



fc2 



-1 

2\k\ 



1 



1 



Po-Po + \k\-ie p'„ - pa - \k\ + : 



(112) 



and closing for each part of ( |112| ) the integration path in order to leave the 
photon's pole outside, with the result 



Hp' J) 



2iTT 



1 



1 



E2-P0 + \k\ - te El 



Po + \k\ 



(113) 

which is the expression obtained in time-ordered perturbation theory |34], and 
again a rearrangement of (109). 

If we include now the four- vertex graphs, we must compute 



V 



(114) 



of which we already have < > . If we choose to split G° into , we 

get, keeping only the terms with zero or one 



V : 



< K^G^K'^ > 



5S 



+ [1- < .gO] < + <K^ [1 - 5° < ] • (115) 

The second line will not contribute to the first-order energy shift of a perturba- 
tion calculation starting with < . If we start with a Coulomb potential, 
then < will provide perturbations in (pc) and the second line per- 
turbations in (pc) . The two last terms of the first line are in (pc) , but 
their higher-order contributions cancel mutually, leaving contributions in 
(pc). 



4.2 Three fermions. 

For the three-fermion problem, we can write 

(l) = g°V(j)^ V ^< K > + < KG"K > - < K > g" < K > 
with the definitions 

<A> = « G"AG° » 



(116) 
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(117) 

5" = ^ , G"^G°G°Gl (118) 

Po- Ei^ E2- Ea + ie' 123 \ ) 

K = K\2 {Gtr' + K23 (G?)-i + A'31 (G°)-i + Ki23 (119) 

dpo = S{Po - pio - P20 - P3o) dpio dp2a dpao- (120) 

We introduced the notation <<>> in order to exhibit the operators 

below. The Kij wih be given by (107). We shaU here neglect the three- 



fermion irreducible kernel which begins with 6-vertex interactions only. The 
two- and 4-vertex diagrams contributing to V are given in figure 2 (we draw 
only L^jCs^B^ and L^i, as the other ones can be obtained by two-fermion 
permutations and/or three- fermion circular permutations). The computation 
of these diagrams is tedious but straightforward. Each external fermion line 
represents a propagator {piQ—Ei + ie)^^ and each internal fermion line a prop- 
agator {piQ — hi + iehi)^^. We could integrate clockwise with respect to the 
energies on the lines 1 and 3 (for the diagrams drawn in the figure), replacing 
P20 by Po~Pio~P30- The leading contributions will come from the four residues 
of the positive-energy fermion poles on lines 1 and 3, while higher-order contri- 
butions will appear when one or several of these poles are replaced by a photon 
pole or by an internal fermion 2 pole (in G3 the residues on the internal 1 and 
2 lines are both leading terms, if we close the integral on piQ clockwise, but 
they cancel mutually at leading order). 

For the sum of the (12)+3 unconnected graphs, the residue of the pole in 
(p30— i?3 + ie)^^ will give the (12) potential, with the (12) energy P120 replaced 
by Pq — E3 (if we close the pso integration path clockwise) . The sum of the 
residues of the propagator of the spectator fermion in the unconnected graphs 
will thus give the sum of the three two-fermion potentials: 

K = < K'[2 >Pi2o=Po~E3 + < K23 >P23o=Po-Ei +< -ft^Ji >P3io=Po--B2 ■ (121) 

In graphs like G3 and B3, however, the internal fermion propagators will also 
provide poles inside the pao integration path, unless we neglect the negative 
energy parts. These contributions (which can be neglected at order ) appar- 
ently violate cluster separability, even if we do not choose the reference frame as 
being the global rest frame and left it unspecified. They do not vanish indeed 
when all interactions between particle 3 and the (12) cluster are "switched off'. 
However, such a "switching off' is of course not compatible with the existence 
of a global bound state, and we have shown above that our 3D equation gives 
anyway the correct cluster-particle scattering amplitudes. 

The connected graph L13 is more specific of a three-body problem. It is 
given by 

ii3 = < G° i^f2 > - < K^, (G?)-i > 5" < K^^2 (Glr' > 
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j^O x^O /^O T^L y^O T^L ^ 

— ^< l^i iV23 ^2 ^12 *-'2 ^-'S ^ ^ 

- « G? G° G° i^2^3 G^ G° » g" « G° G° iff^ G° G° Gg » . (122) 
Each Gi and G3 can be splitted into 

G° = Gf + Gf , Gf = -2nrA+SipM~E,)P,. (123) 

Let us consider the initial G§ and the final G" in ii3. There is a mutual 
cancellation between the two terms of L13 when these propagators are both 
replaced by their S part (this suppression of the leading contribution is indeed 
at the basis of our reduction method). If we keep only the one-G^ contributions 
to L13, we get 

^ ^ r^^ /^o Ti^^ r^^ /^o r^R \ \ 

~ << U2 i\23 "^2 ^12 "-^1 "^2 "^3 -^-^ 

^ ^ r^R /^O i^L r^O T^L y^(5 y^O r^b ^ ^ 
<<. •-^2 "-^3 -'*-23 "^2 ^12 "-^1 '-^2 '-^3 ^ 

< >n + < ^2^3 >3? 5° < ^^2 >ii (124) 
where the indexes 11 or 33 indicates which G" is replaced by G\ or Gf . 
The two first terms of ( |124| ) will contribute in (pc). The two other terms 
will contribute in o? (pc) and the facto rs inside them will be given, up to 
permutations, by one of the three terms of (108- 109| ), with the spectator fermion 
on its positive-energy mass shell. As in the two-fcrmion problem, however, the 
contributions of these terms will be of higher-order in a perturbation calculation: 
changes of indexes < >ii —> < K^j >jj correspond to changes in (pc) in 
(124). If we isolate a leading contribution in (pc) we can built combinations 
like 



L 



13 



^ -"^23 ^33 



[1 - (< K^:, + < K^, + < K^^ 5" ] <K 



12 ^ 



(125) 



coming from L^^B^^Li^ and L23. 

If we perform the energy integrations of ( 124 ), we found, for the first and 
the last terms 



L 



(1+4) 



13 



{P'P] 



(2 7r)e 



A+(g2)(l-al•a2)A^ 
1 



(126) 



1 



2\kl 



Pn-E[-E2m-E'^-\' 
P0 — E1 — E2 — E3 



{E[-Ei)^-\ku\ 



E3-E3- 



Pa — Ei — E2 — E'n—\k'2 



+ 1 



231 



(127) 



where we neglected the A2 (52) part of the internal G2(<72) propagator and its 

['2 13"] 

residue. The part of L13 in obtained by the permutations 



I k 



|fcl2 



^1 



E'l <-> £'31 E2 



E'2. 
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4.3 Gross' spectator model. 

When the three-fermion irreducible Bethe-Salpeter kernel is neglected, the global 
propagator becomes a sequence of two-fermion interactions, for which the third 
fermion is a spectator. In Gross' spectator model |2l], the spectator fermion 
is put on its positive-energy mass shell (we could say that we keep only the 
Gf part of the corresponding propagator G"). When a (12) interaction, for 
example, is followed by a (23) interaction (first part of figure 3), the fermions 1 
and 3 between them are put on their positive-energy mass shell. The neglected 
Gf parts and the neglected three-body irreducible parts of the kernel are then 
treated as higher-order corrections. The result is a set of 3D equations for par- 
tial transition operators, which can not be reduced to a single equation as the 
initial and final "spectator" fermions are not the same in all terms. 

In the integrating propagator based reduction, all terms of the transition 
operator are treated in the same way. We expect of course that the most impor- 
tant contributions to the integrals in the relative energies will be the residues 
of these propagator's poles which put the spectator fermions on their positive- 
energy mass shell. The 3D potential V differs from the integrated transition 
operator < T > by the presence of the counter-terms in — > g'^< , as shown in 
figure 3. If we put the spectator fermions on their positive-energy mass shell in 
these terms too, then the two diagrams of figure 3 cancel mutually, and, finally, 
there remains only the unconnected contributions ( |l2l| ) to the potential. This 
is in fact the kind of potential we proposed in ref. |19|| (in which there is also 
a comparison with Gross' spectator model). The Gf parts of the free propa- 
gators of the spectator fermions and the irreducible three-body terms provide 
corrections to it. 

4.4 Coulomb gauge. 

In Coulomb's gauge, we would make the replacement 

^-/3i/32(7r72) -> ^ - TT^("i-"2)^ (128) 



id,-d,f = {d,-a,) - (""'l^""') . (129) 

k 2 

In (128), the relative-energy dependence lies in the second (transverse) term. 
The contribution of this term to the energy begins at a^, as it connects the 
positive-energy components with the negative-energy components. At this or- 
der, we can replace fc^ by — fc^. In the two-fermion problem, we get 

Iip\p) = (130) 
V{p',p) = -1 ^A++(p')(l-(«i-«2)^)A++(p). (131) 
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In the three- fermion problem, up to order a^, we smiply must add the three 
Coulomb and the three Breit potentials: 

+ V23 ip'23,P23 ) S^'ip'.-pi) + 1/31 (P3l,P31 ) S^p'^-P2)- (132) 

The contributions which are specific to the three-body problem will thus appear 
beyond a"*. 

5 Conclusions 

Our integrating propagator-based reduction of the Bethe-Salpeter equation pro- 
vides a quite straightforward way of writing the different terms of the 3D po- 
tential for any N-particle bound state equation. The continuum-dissolution 
problem is automatically avoided. Our reduction provides a link between the 
kernel-based reductions (especially that of Phillips and Wallace [^) and the 
equal-time retarded full propagator method of Logunov and Tavkhelidze [ p2[ 
and Kvinikhidze and Stoyanov j2^, Although first developed for the two- 
fermion bound state problem, it has been naturally generalized to the N-particle 
bound state problem and is also valid for the 2-body and particle-cluster scat- 
tering problems (even if one may prefer another approach in this case). The 3D 
equation provides a (in principle) indefinitely improvable approach of the exact 
values of the measurable quantities, despite the fact that the individual terms 
of the series giving the 3D potential are neither Lorentz covariant nor all cluster 
separable. 
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Figure 2. The lowest-order three-fermion graphs. 




Fieure 3. Gross' snectator model. 



